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Why people use Level 3 products? 

• Satellite Level 2 data are difficult to work with 
because of: 

- Complex formats 

- Complicated projection (swath) 

- Data volume 

- Number of files 

• Level 3 products are widely used by modelers, 
application users, climate change scientists 

• Level 3 data are easy to use ... but also are easy to 
misuse 
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Why something might go wrong with 

Level 3 products? 

• Usually, Science Teams are tasked to produce & validate Level 2 data 

• Level 3 products are treated mostly as just imagery, to assess gross 
features and variability of geophysical parameters 

• Usability of L3 data usually is not a high priority 

• L3 data are constructed differently for different instruments 

• L2 errors usually are not propagated to L3 

• At best, standard deviations (mostly reflecting variability within a grid 
box), sometimes pixel counts and quality histograms are provided 

• The L3 "validation", in most cases, is done by consistency checking 
and comparing with L3 data from other sensors or models 

• No consistent efforts to characterize & quantify L3 uncertainties 
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So, what do we need to do? 

Despite Science Teoms/Dota Providers not encouraging 
using 13, these products are widely used, so their quality 
and differences between them need to be addressed 

• Instrumental issues: measurement precision, differences in 
calibration and instrument sensitivity, changes or drift in 
calibration or sensitivity over time, etc. 

• Fundamental differences in the retrieval algorithm method 
(multispectral vs. multi-directional) and assumptions: aerosol 
models employed, wavelengths used, cloud filtering, surface 
reflection handling 

• Observational issues: rapidly varying cloud cover, viewing 
angles and conditions, time(s) of observation, changing 
surface characteristics, etc. 
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Differences in L3 from different sensors, cont. 

When comparing data from different sensors, it is important to 
understand and (where possible) consistently process the data. 

L3 data processing: 

• Spatial and temporal binning (L2->L3 daily): 

- Measurements (L2 pixels) from one or more orbits can go into a 
single grid cell different within-grid variability 

- Different weighting: pixel counts, quality 

- Thresholds used, i.e., > 5 pixels 

• Data aggregation (L3D L3monthly regional global): 

- Weighting by pixel counts or quality 

- Thresholds used, i.e., > 2 days 

While these algorithms have been documented in ATBD, reports and 
papers, the typical data user is not immediately aware of how a given 

portion of the data has been processed. 
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Level 3 dataday definition 


Here we address one more aspect of potential 
incompatibility between Level 3 data coming from 
different sensors: the dataday issue. 


Let's compare Level 3 daily data between MODIS- 
Terra and MODIS-Aqua: 


8/25/2010 


Gregory Leptoukh 



MODIS vs. MODIS 



Cornel a tion(A£B) (0lJan2003 - 31Dec200S) 

A : M0D03_D5,005 Aerosol Opticol Depth ot 550 nm (unitlees) 
B: MYDOS D3.051 Aerosol Oob'cal Death ot 550 nm 


-0.B -0.4 — 0,2 


MODIS-Terra vs. MODIS-Aqua: Map of AOD temporal correlation, 2008 
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kfYDQR P5-Q51 Aerosol Optical Depth at 
550 nm (un'rtleas) 



AOD MODIS Terra vs. Aqua in Pacific 

Over the dateline Away from the dateline 




R 2 = 0.45 R 2 - 0.72 

RMS = 0.05 RMS = 0036 
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AOD Aqua MODIS vs MISR correlation map 



Cornel a tian(AAB) (01Jan2003 - 31Dec2008) 

A: MYDG8_D3,051 Aerosol Optical Depth ot 550 nm (unitless) 
MIL3QAEl<irc.Q04 Aerosol Optical Depth at 555 nm (Green Bond 
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AOD Aqua MODIS vs MISR correlation map for 2008 

Note: It is basically the same sensor but on different platform, 

but correlation is not that good 
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MODIS vs. MISR on Terra 
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Cornel Qtion(AAB) (Q1Jon2003 - 31Dec20Q3) 

A; MOD03_D3,005 Aerosol Optical Depth ot 550 nm (unitleaa) 
B: MIL3DAEIarc.QQ4 Aerosol Optical Death at 555 nm ( Green Bond 
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MODIS-Terra vs. MISR-Terra: Map of AOD temporal correlation 

Note: a very good correlation globally besides a narrow area 
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MODIS Atmos. Data day definition 


Level 3 daily products are generated by binning Level 2 data belonging 
to one day onto a certain spatial grid according to a dataday definition. 

Dataday might be different for different sensors and sometimes even for the same 

sensor but defined by different teams. 


MODIS Atmospheric products (from MODIS L3 ATBD): 

The Daily L3 product contains statistics computed from a set of L2 MODIS 
granules (HDF files) that span a 24-hour (00:00:00 to 23:59:59 UTC) 
interval. In the case where a L2 parameter is only computed during the 
daytime, then only daytime files are read to compute the L3 statistics. 
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Orbit Time Difference for Terra and Aqua 2009-01 -06 


Terra 


Terra 


TERRA PREDICTED PASSES JANUAR' 


Orbit track from: http://www.ssec.wisc.edu/datacenter 
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Latitude 


*% Max Time diff. for Terra (calendar day) 


Maximum Local Time Difference for Calendar day definition 
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Data day definitions 


OK, now we know where the artifact has come from. Is there a 
way to fix it, i.e., the right dataday definition? 


1. Calendar dataday: all granules between 00:00 - 24:00 UTC: 
MODIS Atmospheric products, OMI L2G - problematic 


1. Local time (pixel-based): uses local date/time for each pixel 
and ensures spatial continuity - good: 

TOMS, AVHRR, AIRS, OMI, MODIS Ocean, SeaWiFS, MERIS 
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Spatial (local time) Data Day definition 


• Each data set contains information for 24 hours of local time, 
e.g., 13:30 

• The gridding starts at the dateline and progresses westward, 
as does the satellite 

• Parts of scan lines that cross the dateline are included in the 
current date data set or the next, depending on which day is 
at the local time/day at that longitude. 

• For Aqua, the p.m. orbit starts at roughly 1:30 Z on the day 
and ends on roughly 1:30 Z of the following day. 


8/25/2010 


Gregory Leptoukh 


Latitude 


Max time diff. between Terra and Aqua 



Longitude 


Longitude 


Calendar UTC (MODIS) dataday Spatial (local time) dataday 


The artifact around the dateline disappears. In other areas, results are exactly the 
same for the (-7, 18) latitude belt. 

At higher latitudes, the additional restriction for one orbit time around the local time 
produces different results for two dataday definitions. 
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Latitude 


*% Removing the artifact in 16-day AOD correlation 



Calendar dataday 
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MODIS Cloud Top Pressure 



Correia tionlA&B) (0lJan2003 - 1 <5 Jan 200 8) 

A: M0D0a_D3*005 Cloud Top Pressure (Day only) (hPo 


B: MYDOS 03*051 Claud Too Pressure (Dov onl 




MODIS-Terra vs. MODIS-Aqua: Map of CTP temporal correlation, Jan 1-16, 2008 
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($& MODIS Terra & Aqua vs. AIRS Cloud Top Pressure 



10Jan2OO0) 


A: M0D0a_D3.005 O ud Top Pr 'jra {Day only) (hPo) 

6- AIRX3STO.OQ5 Cloud top te iperoture a seen, ^ no {CloudTopTemp A) (Kelvin) 




Correlation' i&B) (0Uon2008 - 6Jan20Q8) 

A: M 0 D 0 a _ D 3.0 5 Cloud Top Pressure Ooy only) 
B: MYD08 D3/ it Cloud Top Pressure CL w only) 


Different 
dataday 
definitions lead 
to different 
artifact patterns 


BON 


■ 3 GN 


MS 


60 S 


■&QN 


Correlation (A 
A: AIRX3STD.Q05 Cloud top 
B: MYDOS 03.051 


Correlation 
maps for Jan 
1 -16, 2008 


- 1.0 — 0.6 


M0DIS Aqua us, M JDIS Terra 

-i n n a M ■ no n n M n a i n 

f A 


8/25/2010 


22 


% Conclusions and recommendations 

• The calendar UTC 00-24 (MODIS) dataday definition 
leads to artifacts around the dateline due to At 
between measurements reaching up to 23 hours 

• Spatial (local-time-pixel-based) dataday definition 
insures consistently small At between measurements 
from different satellites, thus removing artifacts 

• Different "packaging" of L2 into L3 makes a difference 

• Dataday effect is just one of many sources causing 
differences between Level 3 products from different 
sensors that need to be assessed and characterized 
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